Origin of the different conductive behavior in pentavalent-ion-doped anatase and 

rutile Ti0 2 
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The electronic properties of pentavalent-ion (Nb 5+ , Ta 5+ , and I 5+ ) doped anatase and rutile Ti0 2 
are studied using spin-polarized GGA+ U calculations. Our calculated results indicate that these 
two phases of Ti02 exhibit different conductive behavior upon doping. For doped anatase Ti02, 
some up-spin-polarized Ti 3d states lie near the conduction band bottom and cross the Fermi level, 
showing an n-type half-metallic character. For doped rutile T1O2, the Fermi level is pinned between 
two up-spin-polarized Ti 3d gap states, showing an insulating character. These results can account 
well for the experimental different electronic transport properties in Nb (Ta)-doped anatase and 
rutile TiQ 2 . 
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Transparent conducting oxides (TCOs) have many ap- 
plications in optoelectronic devices such as flat panel dis- 
plays, organic light-emitting diodes and solar cells. As 
one of the TCOs, Sn-doped ln 2 3 is widely used because 
of its excellent optical transparency and electrical trans- 
port property. [l| However, owning to the high cost of in- 
dium and the increasing demand for high-performance 
TOCs, many efforts have been made to develop new 
TOCs materials. 0, S] Recently, as one potential candi- 
date of TCOs, Nb (Ta)-doped anatase Ti02 has attracted 
lots of attention because of its hig h electrical conductiv- 
ity and optical transparency. [4-11 However, the origin of 



its high conductivity is still controversial. Wan et aZ. |12| 
found that Nb-doped Ti0 2 grown on (0001) A1 2 3 sub- 
strate shows a much larger resistivity than that grown on 
(100) SrTiOs substrate, and hence the Nb diffusion into 
the SrTiOs substrate is thought to lead to the high con- 
ductivity of Nb-doped Ti0 2 . Meanwhile, Furubayashi 
et aZ. |13(| confirmed that Nb-doped Ti0 2 forms a ru- 
tile phase on the (0001) Al 2 03 substrate but an anatase 
phase on the other substrates, and hence they suggested 
that the higher resistivity of Nb-dopcd Ti0 2 grown on 
(0001) A1 2 3 substrates is caused by the formation of 
rutile phase. Interestingly, later experiments further ver- 
ified that Nb (Ta)-doped anatase Ti0 2 is metallic but 
Nb (Ta)-doped rutile Ti0 2 is insulating. 0, El There- 
fore, one may speculate that the Nb (Ta) doping can lead 
to different conductive properties in anatase and rutile 
Ti0 2 , i.e., conductive for anatase phase but insulating 
for rutile phase, [lij 

In principle, a pentavalent dopant, such as Nb 5+ , 
Ta 5+ , and I 5+ , 17|, [l8[ can release one additional elec- 
tron into Ti0 2 than a Ti 4+ , and introduce donor 
levels. [T9L Eoj Although Nb-doped anatase Ti0 2 has 
been studied using standard density functional theory 
(DFT) calculations, [IB-iil the standard DFT calcula- 
tions within either local density approximation (LDA) or 
generalized gradient approximation (GGA) cannot proj 
erly deal with the strong-corrected Ti 3d electrons. [H 
l23l | In the present work, we studied the electronic 



properties of pentavalent-ion (Nb 5+ , Ta 5+ , and I 5+ ) 
doped anatase and rutile Ti0 2 , respectively, using spin- 
polarized GGA+ U calculations which can give a proper 
description of Ti 3d orbitaIs.[l9l [23L Our calcula- 

tions indicate that Nb (Ta, I)-doped anatase Ti0 2 shows 
an n-type half-metallic character, while Nb (Ta, I)-doped 
rutile Ti0 2 shows an insulating character. These results 
give a good explanation for experimentally observed dif- 
ferent conductive behavior in Nb (Ta)-doped Ti0 2 - 

The spin-polarized GGA+ U electronic structure cal- 
culations were carried out using the Vienna ab-inito sim- 

108-atom 3 x 3 x 3 su- 



lfation package (VASP).[25l 126 
percell of the anatase phase and 72-atom 2x2x3 su- 
percell of the rutile phase are used to model Nb (Ta, 
I)-dopcd Ti0 2 , in which a Ti atom is substituted by 
a Nb (Ta, I) atom. The projector augmented wave 
(PAW) potentials are used to treat electron-ion interac- 
tions and generalized gradient approximation parameter- 
ized by Perdew and Wang (PW91) are used for electron 
exchange-correction functional. [23, A cut-off energy 
of 500 eV and a 3 x 3 x 3 fc-point mesh centered at V point 
are used. The lattice parameters and all the atomic posi- 
tions are fully relaxed until all components of the residual 
forces are smaller than 0.01 eV/A. In our GGK+U cal- 
culations, the on-site effective U parameter (U e ff=U- 
J=5.8 eV) proposed by Dudarev et al, is adopted for Ti 
3d electrons. I29i |30I 



To examine the substitutional Nb (Ta, I) doping ef- 
fects on the electronic property of Ti0 2 , we calculated 
the total density of states (TDOS) and partial density of 
states (PDOS) for Nb (Ta)-doped anatase (see Fig[T]) and 
rutile (see Fig©, and the TDOS and PDOS for I-doped 
anatase (see Fig [3]) and rutile (see Fig SJ) . For Nb (Ta)- 
doped anatase Ti0 2 , the calculated TDOS shows that it 
is spin-polarized, and some up-spin-polarized gap states 
extend from the conduction band (CB) into the band gap. 
These up-spin-polarized gap states are located just below 
the CB bottom, and cross the Fermi level, indicating an 
n-type half-metallic character. This is in good agree- 
ment with the experimentally observed excellent conduc- 
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FIG. 1. (Color online) Calculated (a) TDOS and PDOS plots 
for (b) Ti 3d and (c) Nb Ad (Ta bd ) states for Nb (Ta)-doped 
anatase Ti02- The Fermi level is indicated by the vertical dot 
line at eV. 
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FIG. 2. (Color online) Calculated (a) TDOS and PDOS plots 
for (b) Ti 3d and (c) Nb Ad (Ta 5d ) states for Nb (Ta)-doped 
rutile Ti02. The Fermi level is indicated by the vertical dot 
line at eV. 



tive property in Nb (Ta)-doped anatase TiO9.fi- 11 1 The 
PDOS shows the Nb 4d (Ta 5d ) states mix with the Ti 
3d states in the whole CB, indicating that the Nb (Ta) 
dopant forms a strong Nb (Ta)-O bond. It is also noted 
that the up-spin Ti 3d orbitals strongly hybridize with 
the Nb Ad (Ta 5d ) orbitals, and as will be further dis- 
cussed below, the Ti 3d orbitals mostly contribute to the 
up-spin gap states. For Nb (Ta)-doped rutile phase (see 
Fig. [5]), as in the case of doped anatase, Nb Ad (Ta 5d 
) states spread over the whole CB, and two fully spin- 
polarized gap states with a gap about 1.1 eV lie in the 
band gap. Its PDOS shows that these two gap states are 
also mainly contributed by the Ti 3d states, however, 
the Fermi level lies between the two gap states, show- 
ing an insulating character. These calculated results can 
account for the experimentally observed much higher re- 
sistivity in Nb (Ta)-dopcd rutile phase than the anatase 
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FIG. 3. (Color online) Calculated (a) TDOS and PDOS plots 
for (b) Ti 3d and (c) I 5s /bp states for I-doped anatase Ti02- 
The Fermi level is indicated by the vertical dot line at eV. 
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FIG. 4. (Color online) Calculated (a) TDOS and PDOS plots 
for (b) Ti 3d and (c) I bs/bp states for I-doped rutile Ti02. 
The Fermi level is indicated by the vertical dot line at eV. 

In addition, it is worth mentioning that our calculated 
results are different from previous GGA+ U calculations 
done by Morgan el al., 3l| in which the Nb (Ta)-doped 
anatase and rutile Ti02 both show insulating proper- 
ties. This discrepancy may be attributed to the follow- 
ing two reasons: (1) Morgan et al. applied a U parame- 
ter of 4.2 eV for Ti 3d electrons, which is much smaller 
than the used value of 5.8 eV in this work. (2) In our 
GGA+ U calculations, all the degrees of the freedom for 
Nb (Ta)-doped Ti02, including lattice parameters and 
all the atomic positions, are fully relaxed. In contrast, in 
Morgan et al.'s calculations, only the internal degrees of 
freedom are allowed to relax. Therefore, the difference 
in U values and structural optimization methods may be 
responsible for the different electronic properties of Nb 
(Ta)-doped Ti02 in our calculations from those of Mar- 
gan et al. In fact, Orita also found a metallic character in 
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FIG. 5. (Color online) Calculated spin density distribution of I-doped anatase and rutile Ti02 



Nb-doped anatase TiC>2 using GGA+U calculations. [32j 
in which all the atomic coordinates and lattice constants 
are optimized. In particular, recent hybrid density func- 
tional calculations also show that Nb (Ta)-doped anatase 
Ti0 2 is metallic, while Nb (Ta)-doped rutile Ti0 2 is 
semiconducting. |33J In summary, these related studies 
further confirm our GGA+ U calculations. As a conse- 
quence, we can conclude that the conducting character 
of Nb (Ta)-doped anatase and the insulating character of 
Nb (Ta)-doped rutile are their intrinsic properties. 

For I-doped anatase and rutile Ti02, similar to the case 
of Nb (Ta) doping, an n-type half-metallic character oc- 
curs in I-doped anatase phase but an insulating character 
occurs in I-doped rutile phase. The calculated TDOS and 
PDOS are shown in Fig [3] for anatase phase and Fig [4] 
for rutile phase. However, different from that of Nb (Ta) 
doping, a double filled gap state mostly consisting of I 
5s orbital appears in the band gap, which is located just 
above the valence band maximum. This indicates that I 
dopant exists as I 5+ (5s 2 5p°) in TiC>2, which is consistent 
with the standard GGA calculations. fl7j 

To understand the origin of the different conductive 
behavior associated with the spin-polarized Ti 3d states 
in Nb (Ta, I)-doped anatase and rutile Ti02, we take I- 
doped Ti02 as an example to show its three-dimensional 
spin density distribution in Fig. El For I-doped anatase 
Ti02, its spin density mostly comes from the four equiv- 
alent second-nearest Ti ions around the I dopant (see in 
Fig. [5^). These four equivalent Ti ions share one electron 
donated by one 

J5+ 

ion with their Ti 3d orbitals, and 
thus produce a total spin magnetic moment of 1.0 \xb- 
Therefore, to a first approximation, these four equivalent 
Ti ions should exist as Ti+ 3 - 75 (d 1 / 4 ). Actually, exper- 
imental core-level photoemission spectra measurements 
showed a minor peak of the binding energy below that 
of Ti 4+ ion, [2(| and this chemical shift corresponds to an 
increase of the valence electron density on Ti 3 ti orbitals. 
For I-doped rutile Ti02, in contrast, its spin density is 



mainly contributed by the two equivalent second-nearest 
Ti ions (see Fig[5jD). These two Ti ions share one elec- 
tron and produce a total spin magnetic moment of 1.0 
fig , and thus one can assume that the two equivalent Ti 
ions exist as Ti +3 5 (d 1 / 2 ). As a result, it is expected 
that a lower binding energy of Ti +3 - 5 ions than that of 
rpj+3.75 j ong can kg observed in Nb (Ta, I)-doped ru- 
tile Ti02 through experimental core-level photoemission 
spectra measurements. Furthermore, the increasing of 
the electron density on Ti 3 c? orbitals directly leads to the 
splitting between Ti 3d occupied states and unoccupied 
states, which is responsible for the insulating character 
of I-doped rutile Ti02- Similar spin density distributions 
also occur in Nb (Ta)-doped anatase and rutile Ti02- 

In summary, we studied the electronic properties of Nb 
(Ta, I)-doped anatase and rutile Ti02 by spin-polarized 
GGA+ U calculations. In doped anatase Ti02, the Fermi 
level is pinned in some up-spin-polarized gap states near 
the CBM, showing an n-type half-metallic conductive 
property. In doped rutile, in contrast, two localized states 
with a gap about 1.1 cV are introduced in the band gap, 
and the Fermi level lies between them, showing an insu- 
lating character. Therefore, to prepare the Ti02-based 
TCOs through pentavalent-ion-doping, it is essential to 
avoid the phase transition from anatase to rutile. Our 
theoretical calculations may provide some useful guid- 
ance to develop Ti02-based TCOs. 
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